Decay of a plasmon into neutral modes in a carbon nanotube 
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We evaluate the rate of energy loss of a plasmon in a disorder-free carbon nanotube. The plasmon 
decays into neutral bosonic excitations of the electron liquid. The process is mediated either by 
phonon-assisted backscattering of a single electron or Umklapp backscattering of two electrons. To 
lowest order in the backscattering interactions the partial decay rates are additive. At zero doping 
the corresponding decay rates scale as power-laws of the temperature with positive and negative 
exponents for the two mechanisms, respectively. The precise values of the exponents depend on the 
Luttinger liquid parameter. At finite doping the decay rates are described by universal crossover 
functions of frequency and chemical potential measured in units of temperature. In the evaluation 
of the plasmon decay, we concentrate on a finite-length geometry allowing excitation of plasma 
resonances. 



I. INTRODUCTION 

Electronic properties of single wall carbon nanotubes 
(CNTs) have been studied intensively in the DC limit by 
means of electron transport^, by optical excitationg^HS!, 
and by angle-resolved electron energy loss spectroscopy 
in the UV region^. CNT properties in the intermediate 
(THz) band of frequencies are harder to access experi- 
mentally and consequently they attracted less attention 
in theory. The dominant excitations in that frequency 
range are plasmons. Initial experiments probing them 
by means of time-resolved measurements were reported 
recentlji^. Here we are interested in the intrinsic mech- 
anisms of the decay of plasmons into other excitations 
within a metallic (armchair) CNT. 

The existing tunneling experiments provide support for 
the L uttin ger liquid description of the charge carriers in 
CNT^SHSl^ In that description, the collective excitation, 
plasmon, is just one of the eigenmodes with linear spec- 
trum. As such, it does not decay. 

In the conventional Fermi liquid description applicable 
to higher-dimensional conductors, separation of excita- 
tions into plasmons and Fermi-quasiparticles is not exact, 
leading to a plasmon decay by exciting quasiparticle^^^'. 
In contrast, in one-dimensional conductors quasiparticles 
are ill-defined, and plasmon decays into other collective 
excitations of the electron liquid. 

In armchair nanotubes in addition to the plasmon 
mode the spectrum of bosonic low energy excitations of 
electron liquid contains three neutral modes propagat- 
ing with the Fermi velocity. Conservation of energy and 
momentum allows a plasmon decay into neutral modes, 
as illustrated in Fig. [l] However absence of coupling 
between the bosonic modes in the Luttinger liquid de- 
scription precludes the possibility of plasmon decay. It 
arises from backscattering interactions and thus requires 




FIG. 1: (Color online) Schematic picture of a possible decay 
process of a plasmon excitation with momentum q and energy 
e = Uc+q into two neutral mode excitations with momenta 
gi and 52 and energies ei = vpqi and £2 = —VFq2- The 
line with the steep slope represents the plasmon mode with 
velocity Uc+ and the two lines with shallower slopes represent 
neutral modes with velocities ±«f. Conservation of energy 
and momentum requires e = ei + 62 and q = qi + q2- 



consideration of corrections to the Luttinger liquid ap- 
proximation. 

Another qualitative difference between plasmon decay 
in nanotubes from that in higher dimensions is that in 
one dimension the excitations of neutral modes generated 
during the plasmon decay can propagate only in two dis- 
tinct directions. The plasmon decay products propagat- 
ing along the same direction may recombine into a plas- 
mon. Therefore in long devices perturbative treatment 
of plasmon decay, strictly speaking, is not justified. 

In this paper we consider plasmon decay in a finite 
geometry illustrated in Fig. [2j In this experimentally 
motivatecP2l setup a carbon nanotube connected to two 
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FIG. 2: (Color online) Schematic representation of the device. 
A carbon nanotube (CNT) connected to two metal leads is 
subjected to ac-electric field. 

metal leads is subjected to THz electric field. The finite 
length of the device enables us to account for the plas- 
mon decay processes using perturbation theory in the 
backscattering interactions. We investigate how the de- 
cay processes affect the plasmon resonances and modify 
the device conductance. 

The paper is organized as follows. In Sec.|ll]we present 
the Luttinger liquid description of electron liquid in arm- 
chair carbon nanotubes. In Sec. IIIII we obtain the ac 
conductance of the device in the absence of plasmon de- 
cay into the neutral modes. In Sec. |IV| we discuss the 
physical mechanisms of plasmon decay and study their 
influence on energy dissipation in the ac conductance of 
tlie device. Our results are summarized in Sec. [V] 

II. LUTTINGER LIQUID DESCRIPTION OF 
ARMCHAIR CNTS 

Near the Fermi level the electron band structure 
in armchair nanotubes consists of two spin degenerate 
bands crossing at two points, k = iQ, in the Brillouin 
zone. Near the crossing points the spectrum has linear 
Dirac-likc dispersion, see Fig. [Sj For undoped tubes the 
tlic Fermi level goes through the crossing points. Setting 
h — 1 we write the non-interacting electron Hamiltonian 
near the Fermi level as 

Hq = -ivp^ / dx^+^^{x)rdx^ara{x)- (1) 

ara 

Here Wi? w 8 x 10^ ms~^ is the Fermi velocity, a = ±1 is 
the valley index, r = ±1 represents right- and left- mov- 
ing electrons, and a is the electron spin. The electron 
creation/annihilation operators are denoted by tp'^ /ij}. It 
is important to note that armchair CNT are symmet- 
ric about a plane containing the tube axis and bisecting 
the AB bonds that point along the circumference. As 



a result the electron wave functions can be classified by 
parity P = —ar = ±1, as indicated in Fig. |3] The wave 
functions with positive/negative parity are respectively 
even/odd in the AB sublattices. 

We consider an {N, N) armchair CNT, which has N 
unit cells around the circumference of the tube, and use 
1/N as a small parameter in the theory. For most syn- 
thetic CNTs > 8, and this is a good approximation. 

For <C 1 the forward scattering part of the e- 

e interaction is much stronger than the backscattering 
part, the latter being small by 1/JSl^. The former can 
be written as 

= y / dxn'^ix), n{x) = ^'il}+^^{x)ipara{x). (2) 

ara 

We assume that the Coulomb interaction is screened, 
e.g. by a gate at a distance d, which is assumed to 
be short compared to the length of the tube, L but 
longer than the CNT radius R. The latter can be ex- 
pressed in terms of the graphene lattice constant a as 
R — \/3Na/2TT. Elementary electrostatics leads to the 
estimate gp « 2e^ ln{d/R). 

Thus to zeroth order in 1/N backscattering is absent, 
an d elec trons in an armchair tube form a Luttinger liq- 
uicP^Esl. The Hamiltonian is bosonized by the standard 
procedure (see, for example, Refs. I17)18p 

Fa a 

^ara{x) = exp{i[ld aa{x) - r ^ aa{x) +rkFx]} . (3) 

Here Faa- are the Klein factors, ^ is a short distance cut- 
off on the order of the radius of the nanotube, hp is the 
Fermi wavevector counted from the Dirac point and $ 
and Q are bosonic fields satisfying the commutation re- 
lation [^aa-ix),Qa'a-'ix')] — —i5aa' Sacr'Q{x — x') , whcrC 

9[x — x') is the step function. 

Introducing the symmetric and antisymmetric combi- 
nations of charge and spin modes in the different valleys, 

^aa = [$c+ + a$c- + 0-$,+ + a(7$,_] /2, (4) 
Qaa = [6^+ + aOc- + 0-6^+ + aaQ,4 /2, 

we rewrite the Luttinger liquid Hamiltonian of the CNT 
including forward scattering as 

H,+Hp = T.^J dx[Kr\d^<^,f + K,{d^e,n (5) 
i 

Here j = c±, s±, and the velocity and Luttinger parame- 
ter of the j-th mode are related via Uj = vp/Kj. All four 
bosonic modes have linear spectrum, Sj — Uj\k\. For the 
three neutral modes j = c— ,s±, the Luttinger param- 
eter Kj is equal to unity and their velocities coincide 
with the Fermi velocity. For the charge mode, j = c+, 
tlie Luttinger parameter is significantly smaller than one, 
Kc+ = + 4,gp/TTVp ~ 0.2. For the reason, the plas- 

mon mode has a much higher velocity than that for the 
three neutral modes. Below we denote the plasmon ve- 
locity by u = vp / Kc+. 
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function by the usual Kubo formula, 




'{[3{x,t),j{x',0)]). (7) 



We assume that the electric field (i) varies smoothly on 
the scale of the lattice constant within the tube and (ii) 
preserves the reflection symmetry of the tube. Because of 
(i) it will not scatter electrons between the valleys, which 
involves momentum transfer ~ 2Q. Due to ii) parity con- 
servation also forbids backscattering of electrons within 
the same valley. Thus the electric field couples only to 
the smooth part of the current density, which can be ex- 
pressed in the bosonized form as j = 2^dt^c+- Therefore 
we can express the conductivity in terms of the Green's 
function of the field. 



FIG. 3: (a) Schematic picture of an armchair tube. Full and 
open circles denote atoms in the A and B sublattices. The 
arrows show atomic displacements in the TA and LO phonon 
modes, (b) Free electron spectrum of armchair carbon nan- 
otubes: + and — denote the parity of the bands. The two 
valleys are labeled by the index a — ±1. 



In the absence of decay processes the energy of a plas- 
mon excited in the device dissipates only due to its emis- 
sion into the leads. The power dissipation due to this 
mechanism is considered in the next section. 



III. PLASMON RESONANCE IN THE 
ABSENCE OF INTRINSIC DISSIPATION 

We consider the setup shown in Fig. [2] The plasmons 
excited in the tube get partially reflected from the junc- 
tions with the leads. The plasmon reflection processes 
can be described by a one dimensional model . Eg. ([s]) 
with spatially dependent Luttinger parameters^fflt^H (de- 
fined as follows. For the charge mode the Luttinger pa- 
rameter is equal to Kc+ inside the tube, < a: < L, and 
unity in the leads. To zeroth order in 1/N the Luttinger 
parameters in the neutral modes are equal to unity both 
in the tube and in the leads. 

We denote the THz electric field by E^{x). Its distri- 
bution along the wire depends on the details of the device 
geometry and is a solution of the proper electrodynamics 
problem which is beyond the scope of this paper; here 
we assume (x) is known. The dissipation rate of the 
electric field energy in the device can be written as 



W{uj)^l^Ref dxi ( dx2E^{xi)E^{x2)<J^{xi,X2), 
^ Jo Jo 

(6) 

where ai^{x^x') is the nonlocal conductivity, which can 
be expressed via the retarded current-current correlation 



a^{x,x') 



4e 



^ G!^{x,x'), 



(8) 



where 

G^{x, x ) 



-I I die'"*([$,+ (x,i),a>,+ (a:',0)]) (9) 
/o 



is the retarded GF of the plasmon mode. 

We denote the Green's function of the plasmon mode in 
the absence of backscattering by Gq^{x,x'). It satisfies 
the differential equation 



dr. 



Vp 



K^{x) 



Vp 



dA+ — \Gl^{x,x')=^8{x-x'), (10) 



with the boundary conditions: (1) Gq^{x, x') is continu- 
ous at X = 0, L and x = x'\ (2) \vp lK'^(x)\dxGQ ^^{x,x') 
is continuous at a; = 0, L. For our model (K{x) equal 
to Kc+ inside the tube and unity outside) the Green's 
function can be expressed in terms of the reflection am- 
plitude, r = i+k"^ ' plasmon from the junction 
with the lead. 

The dissipated power in Eq. ^ depends only on the 
Green's function inside the tube. In this region for < 
X < x' < L it is given by 



GLi^,x') 



ttK, 



c+ 



(x' 



-i — {x' — L) 



(11) 



In the region < x' < x < L the Green's function is 
obtained by exchanging x and x' in the right hand side. 

The Green's function above has poles as a function of 
frequency lo. The real part of the pole location gives 
the resonance frequency a;„ = wku/ and its imaginary 
part gives the resonance half-width, 7^. The plasmon 
reflection amplitude from the junctions with the leads is 
determined by the ratio of the plasmon velocities in the 
tube and in the lead, r = !" , . This gives the escape 



half-rate 7z = f; log jz^^ ■ In the limit of high charge 



4 



stiffness, Kc+ ^ 1, it becomes much smaller than the 
spacing between consecutive resonances and independent 
of the interaction strength, 7^ ~ 2Kc+j^ ~ 2vp/L. 

The dissipated power is given by Eqs. ([s]), and 



(111 and depends not only on the properties of the tube 



but also on the spatial distribution of the electric field (as 
mentioned above, the latter depends on the geometry and 
is assumed known). 



IV. PLASMON DECAY IN THE DEVICE 

The energy dissipation described by the ac conduc- 
tance obtained in Sec. Ill corresponds to the emission 

The THz electric field 
which af- 



of plasmon waves into the leads, 
excites plasmon modes of the electron fluid, 
ter a few reflections from the contacts with the leads are 
radiated out into the leads. Thus the energy absorbed 
from the THz radiation is re-radiated into the leads in 
the form of a plasmon flux. 

In the presence of backscattering interactions, there 
is an additional contribution to the dissipation energy in 
the device. The plasmon waves generated by the THz ra- 
diation can decay into neutral modes before escaping into 
the leads. In this dissipation channel the absorbed energy 
is carried into the leads by the neutral modes. We discuss 
the microscopic interactions responsible for the decay of 
plasmons into the neutral modes in Sec. |IV A| and study 
how these processes modify the device conductance in 
Sec. UVBl 



A. backscattering processes leading to plasmon 
decay 

The possible backscattering processes obey certain 
symmetry selection rules. In particular parity conser- 
vation forbids single electron backscattering in armchair 
carbon nanotube due to e-e Coulomb interaction, which 
also attributes to the Klein paradox in graphene^^ but 
two, or more generally, any even number of particles can 
be back-scattered. The electron interaction Hamiltonian 
contains such processes but their dimensionless coupling 
constants are small in 1 /N and can be treated perturba- 
tively. Having in mind temperatures well below the gap 
between the one-dimensional subbands {vp/R ~ leV) we 
retain only the most relevant backscattering interactions. 
Namely the processes, which sca tter tw o right movers 
into two left movers or vice-verscP^'i^^. Although for 
the Brillouin zone shown in Fig.jsjb) such scattering con- 
serves quasimomentum and thus corresponds to normal 
processes, we refer to it as Umklap p scatter ing, following 
the accepted nanotube terminolog} *^"^ * ^^ ! ^ . 

In the bosonized representation the Umklapp Hamil- 



tonian density can be written a^ 
cos(2$c-i--4fcFa;) 



'Hu — — 



2(7re)2 
(53-5i)cos(2$3+) + 
5i[cos(2$e-)-cos(2$,_)]}. 



{gs cos(2e,_) + 



(12) 



where the coupling constants originate from the short 
range part of the Coulomb interaction. They can be ex- 
pressed in terms of the matrix elements of Coulomb in- 
teraction between electrons residing in the A and B sub- 
lattices, VAA{q) and VAsiq), as follows: 33 = Vaa{'2Q) + 
VabC^Q) and gi = Vaa{0) - VAsiO?^ Here 2Q is the 
momentum difference between the two Dirac points, see 
Fig. [3) and we assumed that the Fermi momentum due 
to finite doping is small, kp <^ Q. Both gi and g^ arise 
from distances on the order of the lattice constant and 
are on the order of e^/N. 

In the presence of phonons single electron backscat- 
tering can occur because of electron-phonon interaction. 
The phonon modes in carbon nanotubes are classified to 
two types by the parity symmetry: i) transverse acoustic 
(TA), longitudinal optical (LO), and radial optical (RO) 
with negative parity, and ii) transverse optical (TO), lon- 
gitudinal acou stic (LA ), and radial acoustic (RA) with 
positive paritjJ^SEfiEH. The couplings of radial phonon 
modes to electrons are weaker in 1/iV^^ than those for 
other phonon modes. Therefore we neglect their effects 
here. For the remaining couplings, one can distinguish 
between inter-valley and intra-valley scattering. The for- 
mer requires a phonon with a momentum on the order 
of the inverse lattice spacing. Such phonons have high 
frequencies ~ 10'^ K and are frozen out at room tem- 
perature. We thus neglect such processes and consider 
only intra-valley backscattering. Such processes change 
the parity of the electron wavefunction and thus require 
phonons with negative parity, i.e., TA or LO phonons. 
Their displacements are illustrated in Fig. [3j Since the 
frequency of the LO phonon is ~ lOOOif it is frozen out 
in the temperature range we are interested in. Thus be- 
low we restrict our treatment to the backscattering of 
electrons from the TA phonons only. 

The Hamiltonian of the TA phonon mode is given by 



Hta — 



dx 



1 

ST 



{ipT{x)f + ST{d^LpT{x)y 



(13) 



Here ^Pt{x) = ^ pstUt{x) is the dimensionless TA 
phonon field, p the linear mass density along the tube 
axis, ut{x) - the atomic displacement of the TA phonon 
mode, and st ^ 1.4 x 10"* m/s is the phonon veloc- 
ity. The deformation of the TA mode couples to the 
time-reversal invariant operator bilinear in the fermions 

M{x) = -iT.ara^T"4'ira(.x)lka-ra, 



Hp 



dxM{x)gTdx^T{x)- 



(14) 



In the bosonized representation the Hamiltonian den- 
sity for the electron- TA phonon coupling can be written 
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ep 



dxVT 



Y[ cos{^c+-^kFx)coa{^c-)sm{^su) 



sin($c+ - 2kFx) sin (<i>c_) cos($si,) 



where the index rj denotes either the Umklapp (77 = u) 
or the electron-phonon (rj = ep) interaction, and the ker- 
nel IC':^{xi,X2) can be expressed in terms of the retarded 
correlation function of the corresponding backscattering 
perturbation as 



(15) 



where gx is the e-ph coupling constant. Within the tight- 
binding model, the coupling constant can be expressed 
in terms of the derivative of the transfer integral J{r) 
with respect to the bond length r, see Ref. [2E- In an 



armchair carbon nanotube, 



dJjr) 
dr 



and 



4 y/pST 

for {N, N) armchair nanotubes, p = ANM/a, where M 
is carbon atom mass and a is the lattice constant of a 
graphene sheet. 



B. intrinsic energy dissipation in the device and 
plasmon decay 

The backscattering interactions lead to the possibility 
of plasmon decay into the neutral modes. Using the fact 
that the backscattering interactions are small in 1 /N we 
evaluate the decay rate and the corresponding correction 
to the device conductance using perturbation theory. We 
work in the regime where the cumulative plasmon decay 
rate in the device is smaller than the rate of elastic escape 
of plasmons into the leads. In this case the correction 
to the device conductance due to plasmon decay can be 
treated perturbatively. 

Both the Umklapp and the electron-phonon backscat- 
tering interactions are expressed in terms of exponentials 
of the charge mode, <&c+- As a result, their contribution 
to the plasmon decay and the device conductance can be 
evaluated along similar lines. 

To second order in perturbation theory there is 
no cross-correlation between the Umklapp and elec- 
tron phonon backscattering processes, and the retarded 
Green's function of the charge mode can be expressed as 

V/ dxi I dx2G^^{x,Xi) X 
]C^ixuX2)G^Jx2,x'), (16) 



K.':;ix,x') = -lal dt[e"^mnr,{x,t),n^{x\0)])o- 



dc{[nr,ix,t),n,,i^,o)])oSix - x') 



(17) 



Here is the coefficient in the exponential dependence, 
~ (e'°''*'=+ -I- /i.e.), of the corresponding backscattering 
perturbation on the charge field $c+- The derivation of 
Eqs. ([l6| an d ^ is outlined in Appendix |Aj The first 
term in Eq. ( 16 1 is the Green's function of the plasmon 
in the Luttingcr liquid approximation, Eq. (11). The 



remaining terms represent correction due to backscatter- 
ing interactions. The first term in the right hand side 
of Eq. (171 corresponds to the usual RPA-hke diagram. 



the second term comes from the tadpole diagram and de- 
scribes elastic scattering. The average (. . .)o is performed 
with respect to the quadratic Luttinger liquid Hamilto- 
nian, Eq. (§. 

In the spatially inhomogeneous setup considered here 
the kernel /C5^(a;i, 2:2) depends not only on the coordi- 
nate difference but also on the center of mass coordinate 
{xi + X2)/2. Its evaluation leads to rather cumbersome 
expressions. The situation simplifies in the physically 
relevant regime, where the characteristic thermal wave- 
length of the neutral modes, Lt = vp/T, is much shorter 
than the size of the device L, while the wavelength of the 
plasmon resonance is on the order L. 

In this case the Green's functions of the plasmon mode 
in the second term of Eq. ( 16 ) depend on xi and X2 on 
the scale of wavelength of the plasmon resonance, ~ L. 
On the other hand, because backscattering operators de- 
pend exponentially on the neutral modes, Hri{x,t) ^ 
exp(ia^$j(a;,t)), their correlator ]C'^{xi,X2) falls off ex- 
ponentially for — X2\ > Lt, which is much shorter 
than the length of the tube L. This is especially appar- 
ent in the case of a uniform wire for which the correlators 
of exponentials of bosonic fields have a simple form^^. 



^g»a$, (.,t)g-»a*, (0,0) ^ n(X, , a; X, t) 



{sinh[7rT(iCja;/uF - sinh[7rT(Xja;/wF + t-)]} 



(18) 



where t_ = t — ie. still decays exponentially for \xi — X2\ > Lt- Therefore 

In the inhomogeneous situation the kernel /C5i'(xi, a;2) tl^e main contribution to the Green's function correction 
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in Eq. ( 16 1 comes from the region — X2 1 < Lt- At such 



separations only the fluctuations of bosonic modes with 
wavelengths less than Lx contribute to the correlator. In 
the regime of interest, Lt L, such short wavelength 
vibrations can be described semiclassically and are insen- 
sitive to the boundary conditions at the junctions, which 
describe reflection processe^'^^. Thus we come to the con- 
clusion that in the bulk of the nanotube, i.e. with the 
exception of the regions of size Lt near the junctions the 
correlator /C^(a;i,a;2) is the same as in a homogeneous 
tube, and can be replaced by the local kernel 



IC':;^{xi,X2) S{xi - X2)/C5^((7 = 0), 



(19) 



where 



/C^(g = 0) = -ia 



00 poo 
2 I I ^Uju.t 



dtie'"^' - 1) 



x([H^(x,t),H„(0,0)])o. 



(20) 



This observation enables us to express the correction 
to the dissipated power due to decay of plasmons into 
neutral modes in terms of intrinsic parameters of uniform 
nanotubes. Accordingly, below all correlation functions 
of backscattering operators will be evaluated for infinite 
homogeneous nanotubes. 



With the aid of Eq. ( 19 ) the power dissipation rate in 
Eq. (|6| can be written in the form 



Wiuj) = Wo{Lo) + lm}_^K:';^{q^O)F{u). (21) 



Here the first term corresponds to the dissipation rate 
due to the escape of plasmon to the leads and is given by 



7r2 , 



dxi 



dxo 



xlmE^{xi)El{x2)Gl^{xi,X2). (22) 



The second term in Eq. (21 ) corresponds to the intrinsic 



dissipation due to the electronic Umklapp processes and 
electron-phonon backscattering. It can be expressed in 
terms of kernel for an infinite tube, Eq. (20), and the 
form-factor Fijjj), 



F{uj) 



dxi I dx2 / dx3E^{xi) E^{x3) 



Jo Jo Jo 

^[xi, X2) Gq (2^2 7^3)7 



(23) 



which is quadratic in the THz electric field and depends, 
through its spatial distribution, on the geometry of the 
device. 

Both the electric field and the unperturbed Green's 
function Gq^ of the Luttinger liquid are temperature in- 
dependent. Therefore Wo(w) and the form-factor F{uj) 
do not depend on temperature. All the temperature de- 
pendence of the dissipation rate W{u}) is thus totally de- 
termined by intrinsic properties of a uniform CNT and 
is described by the kernel IC':^{q = 0). 



We note that the kernel JC':^{q = 0) also determines 
the ac-conductivity cri^{q = 0) of an infinite wir m 
the high frequency or high temperature regime, where 
back-scattering interactions can be treated perturba- 
tively. The latter can be expressed in the form 



a{uj,T)- 



(24) 



Here 7,, is the relaxation rate due to the Umklapp (77 — u) 
or electron-phonon (r/ = ep) back-scattering, which can 
be expressed in terms of the local kernel as 



7»7 = -Kc+Uc+JC';'j{q = 0). 



(25) 



Thus with the aid of Eqs. ml, (23) and (25 1 the tem- 



perature dependent part of the power dissipation in the 
device can be expressed in terms of the intrinsic plas- 
mon decay rate in a uniform CNT, "f,-p which is given by 
Eqs. (l25| and (I20I). 



In the remainder of this section we evaluate the in- 
trinsic plasmon decay rate due to various backscattering 
processes. 



1. Umklapp mechanism of attenuation 

We begin by considering Umklapp-mediated plasmon 
decay rate, 7„. The Umklapp Hamiltonian density in 
Eq. ( [l2| ) contains exponential operators of both <i>- and 
0-fields corresponding to the neutral modes. Because 
the Luttinger parameter K for neutral modes is equal to 
unity, the correlators of exponential operators of and 
0-fields coincide. As a result we obtain 



{Hu{x,t)Hu{0,0))o = 



35? + 2.gi 



16(7r^)4 

U{K,+ ,2;x,t)U{l,2;x,t). (26) 



The analogous correlator corresponding to the opposite 
ordering of backscattering operators can be obtained 
from the identity 

(H„(0,0)H„(a;,t))o = {nu{x,t)nu{0,0))*o. 



Using Eqs. (|25|), (|20|), (|18|), and (|26|) we obtain 



In 



2K^ 



+ igl + 2gl - 2.9133 



4r«)2 



T' T 



(27) 



where FuC^, "^^^ ) is a scafing function depending on 
the ratio w/T and vpkp/T. It is given by the following 
double integral. 
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T' T 



T 



dX 



1 1 cos 



/ Akpvp 



X 



(sinh7r(^X + y_) sinh7r(^X - sinh7r(X + y_) sinh7r(X - 11) 



(28) 



where Y- = Y — ie. press the real part of the scaUng function in the form that 

Using the standard deformation of the integration con- is more convenient for numerical evaluation, 
tour, y_ — > y — 2/2, in the above equation one may ex- 



— — sinh — 
K^oj 2T 



dx 



dy 



cos {yLxj/nT) cos {Akpvpx/irT) 



[cosh {Kc+x + y) cosh {K^+x — y)]^"^ cosh (a; + y) cosh (a; — y) 
I 



(29) 



In the limit of high charge stiffness, Kc+ 0, the above 
expression simplifies to. 



vpki 
T 



cosh ( 



AT 



) cosh 



' ikpVF- 

\ 4T 



(30) 

In the temperature regime of inte rest, T ^ lu u/L, 
the function F„(f , ^^) in Eq. (lisl reaches the dc limit 



F„(0, ^ ) , which is purely real and depends only on 
the doping level. For comparison, in Fig. [4] we plot the 
static limit of the scaling function Fu{0, kpvp /T) eval- 
uated numerically at Kc+ — 0.2 using Eq. ([29]) and the 
analytical expression from Eq. (301 for Kc+ — 0. 



At small doping, the scaling function goes to a con- 
stant and exhibits a power law dependence on the 
temperature, 7„ ~ y2Kc+-i^ strong doping, 2kpVF > 
T the scaling function acquires an activated tempera- 
ture dependence. This results in exponential decay of 
7t, y2/fc+-i exp{—2kpvp/T). The physical reason for 
the exponential temperature decay can be understood 
from simple inspection of energy and momentum con- 
servation, which shows that energetically optimal pro- 
cess involves backscattering of two thermal quasiparti- 
cles exactly at the Dirac point k = Q. At strong doping, 
kpvp ^ T the occupancy of such states is exponentially 
small, ~ exp{—2kpVF/T). 



Equation (27) and the subsequent equations express 



the plasmon decay rate 7„ in terms of the somewhat ar- 
bitrary short distance cutoff, ^, which has the order of 
the radius of carbon nanotube. Alternatively, 7„ can be 
expressed in a renormalized form as a function of an in- 
trinsic cutoff-independent low energy scale. In the case 
of Umklapp interactions the role of such a low energy 
scale is played by the Umklapp gap. Namely, the Umk- 
lapp interaction is a relevant perturbation relative to the 
Luttinger liquid fixed point and drives the system to a 



0.25 I 



0.15 - 



0.10 



FuiO, X) 




FIG. 4: Dependence of the scaling function in the static limit, 
-F„(0,a;) on the doping parameter x — ''K^^ for Kc+ — 0.2 
(solid line) and Kc+ = (dashed line). 



Mott insulating state with the gapP^l, 



5« 

Vp 



Here Qu ~ e'^ /N is the strength of the Umklapp inter- 
action. This gap depends on the radius of the CNT as 
R^\-^+^/(^^^'^<--+)^ and is hard to evaluate theoretically due 
to the ambiguity of the coupling constant, but numerical 
results estimate that it's of the order of 1 — 10 meV for 
an armchair CNT of radius ~ 1 nnPH. For the pertur- 
bation calculation to be applicable, the temperature we 
consider should be above the Umklapp gap. Up to a nu- 
merical coefficient on the order unity the plasmon decay 
rate in Eq. (271 can be expressed in terms of this gap as 



A,, 



T 

a: 



2Ka 



F„, 



u vpkp 
T' T 



(31) 



At zero doping, 7,j ~ T for K, 



c+ 



1, i.e., without 
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considering the Luttinger liquid effect. Yet in the Lut- 
tinger Uquid description with strong Coulomb interac- 
tion, K,+ ^ 0, 7„ - A2/r. 

Note that according to Eqs. (24) and (31 1 the Umk- 
lapp processes yield (t(w — 0,7j c>c at zero 

doping and high temperature, which is consistent with 
the result obtained in a previous worlP^. In the case 
K(.+ — ^ 1, i.e., in the absence of Luttinger liquid effects, 
the decay rate due to Umklapp processes 7„ ~ T, and 
the resistivity has linear dependence on temperature at 
high temperature, which is consistent with the pertur- 
bative calculation of plasmon decay rate from Umklapp 
processes in CNTs^"^ We can see that the inclusion of 
the long range interaction significantly modifies the de- 
pendence of the plasmon decay rate on temperature due 
to Umklapp processes. In the case of strong Coulomb 
interaction, i.e., Kc+ 0, the temperature dependence 
of decay rate is strongly modified and becomes 7„ ~ 
at high temperature. 



2. Phonon- assisted mechanism of attenuation 



potentials becomes effectively short range. In this respect 
the phonon-assisted plasmon decay mechanism is similar 
to that induced by a short range disorder potential. The 
temperature dependence of disorder-induced conductiv- 
ity was obtained in Ref. [M (^dis ^ (y-2-^T^K^+-i)/2 ^ 
where {V"^) is the variance of the random potential. In 
contrast to the temperature-independent variance of the 
disorder potential, the variance of the phonon correla- 
tor scales linearly with the temperature. This gives the 
power law ~ T*^^"*"^=+^/^, which is confirmed by the 
calculations below. 

The phonon-assisted decay rate 7ep can be evalu- 
ated similarly to the case of Umklapp processes. Us- 
ing Eq. ( 18 ) we obtain the following expression for the 



correlation function for the electron-phonon Hamiltonian 
densities given by Eq. (15), 



(Hep(x,t)Hep(0,0)) 



COS {2kpx)Dph{x, t) 



(32) 



We now consider plasmon decay rate due to electron- 
phonon interactions, 7ep. Its physical mechanism can be 
understood as follows. Parity-breaking deformations in 
the TA mode cause single particle backscattering, which 
causes the plasmon decay. Due to the presence of the 
adiabatic parameter sx/vp ^ 2 x 10~^ ^ 1 the lat- 
tice deformations are effectively static. Thus the en- 
ergy of the decaying plasmon is transferred into the neu- 
tral modes of the electron liquid. Therefore this mech- 
anism, which we analyze below, may be characterized 
as phonon-assisted plasmon decay. In two dimensions 
it has been studied in Ref. At finite doping single 
electron backscattering involves phonons with momen- 
tum 2kp. With the exception of very low temperatures, 
T < 2kpST, the occupation numbers of such phonons 
are large, ^ T/2kpST- This corresponds to classical 
thermal fluctuations of the lattice. Since the proba- 
bility of such fluctuations is determined by the elastic 
energy ^ exp(— st / dx[d.j;<j>T{x)]'^ /'2T), the one-electron 
backscattering potential ^ dx(f>TAix) (see Eq. (14)) may 
be viewed as (5-correlated in space. More precisely, 
the correlator of backscattering potential ~ dx<t>TA{x) 
becomes exponentially small at distances greater than 
st/T, as is clear from Eq. (33) below. As the electron 



correlators depend on the spatial coordinates on much 
larger scales ~ Lt ^ st/T the correlator of deformation 



where j = c±, s±, and Dph{x, t) = {\/(Pt{x, t)\/ipT{0, 0)) 
is the phonon correlator given by 



Dph{x,t) - ^ 



2sinh^[7rT(^ -|-i^t_)] 

i/=±l I- ^ St 'J 



(33) 



The correlator with opposite ordering can be obtained 
from Eq. (32) using the identity ('Hep(0, 0)'Hep(x, t)) = 

(Hep(2:,t)Hep(0,0))*. 



Substituting Eqs. (l32|), (|18| and ^ into Eq. ^ and 
using Eq. ( 25 ) we obtain the phonon-assisted rate 7ep . 
While carrying out the spatial integration in Eq. (20) 



we may neglect the position dependence in the electronic 
correlation functions and replace them by n(ii'j,l;a; = 
0, t) (this is justified by the small value of the adiabatic 
parameter st/vp). Doing so we obtain 



>2 """^-F 
7ep = 9t ^ 



ttT^ 

Vp 



sxkp 
T 



(34) 

where cjT = gx / y/srvp is the dimensionless coupling 
constant of e-ph interaction, F^pC^, "^j'^ ) is the scaling 
function for the phonon-assisted decay rate. Introducing 
dimensionless variables X = xT/st and Y — Tt we can 
express it as 



' LU Spkp 

T' T 



T 

W JO 



dY dX (e'T^ - 1) 

Jo J-oo 



COS 



2kpST 
T 



X 



sinh^ {tt{X + Y_)) sinh^(7r(X-y_)) 



sinh (ttY-) sinh(— 7ry_) 



(35) 
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Similarly to the Umklapp case, the real part of the above scaling function can be expressed in a more convenient 
form, 

^Fe, ^) = sinh ^ 1^ dy 1^ dx 



cosh^ (x + y) cosh^ (x — y) 



cos ^ cos ^^^^f?^ 

(c0sh?/)(3+^<^+)/2 



(36) 



Prom the above equations it is clear that the phonon- 
assisted decay rate is practically unaffected by doping 
as long as T > 2kpST- Because of the slow velocity 
St of the phonon mode, with the exception of very low 
temperatures, this condition is satisfied for all practically 
accessible doping levels. In this regime Eq. (36) simplifies 
to 



— smh 2 X 



B 



3+Xc 



iuj 3 + K, 



c+ 



27rr' 



In the limit a; <C T, the scaling function goes to a con- 

stant, Fep(0,0) = 2^^B([3 + if,+]/4, [3 + K,+]/4:)/tt\ 
Similarly to the case of Umklapp processes the phonon 
assisted decay rate can be expressed in terms a renor- 
malized form, which does not involve the short distance 
cutoff. In this case the role of the low energy scale for the 
electron-phonon interaction with the TA mode is played 
by the twist Peierls gajPSl, A^p ~ ^ 9t^^~^''*\ which 
depends on the radius of the CNT as ~ ij-[i+2/(i-ifc+)] 
and in the limiting case Kc+ — is on the order of a 
few tens of Kelvin for a (5, 5) armchair CNT. This gap 
is induced by the Peierls instability of the armchair tube 
with respect to a twist distortior(22H20|_ Up to a numerical 
coefficient on the order unity the phonon-assisted decay 
rate can be expressed in terms of the twist Peierls gap 



T 



ep 



l + -F^c+ 



r' T 



In the dc limit and at zero doping, the Luttinger liq- 
uid effects manifest themselves by changing the power 
dependence of or resistivity from 7ep ^ T for non- 
interacting case'^^'to 7ep ~ T 2^ in a Luttinger liquid 
at temperatures above the Peierls gap. Experiments on 
phonon-induced resistivity to datP^Ml have been con- 
ducted on CNTs placed on substrates rather than freely 
suspended CNTs, and exhibit linear temperature depen- 
dence of resistivity, consistent with the non-interacting 
electron approximation, Kc+ = 1. 



V. SUMMARY OF THE RESULTS 

We studied the intrinsic plasmon decay rates in clean 
armchair nanotubes. Plasmon oscillations decay into 



neutral electron modes. Such processes are facilitated 
by electron-electron and electron-phonon backscattering 
processes. In second order of perturbation theory in the 
backscattering interactions the total decay rate is given 
by the sum of partial decay rates due to each mech- 
anism. For undoped CNT, both electron-electron and 
electron-phonon rates reveal power law dependence on 
the frequency of oscillations or the temperature of the 
system with non-integer exponent, that depends on the 
Luttinger liquid parameter Kc+. 

The plasmon decay processes lead to temperature de- 
pendent corrections to power dissipation in the device 
depicted in Fig. [2j Therefore they can be studied experi- 
mentally in a finite device by measuring the temperature 
dependence of power dissipation. In the regime where the 
temperature length Lt = vp/T is shorter than the plas- 
mon wavelength the temperature-dependent part of the 
dissipated power can be expressed in terms of the instrin- 
sic plasmon decay rate 7 in homogeneous nanotube using 
Eqs. ( 2l]|23 ) and (25). Furthermore, in this regime 7 can 
be related to the intrinsic ac conductivity of carbon nan- 
otubes via Eq. ( 24 ) . Thus intrinsic plasmon decay rate in 



CNT can, in principle, be extracted from measurements 
in finite size devices. 

Among the electron-electron interactions resulting in 
plasmon decay the most relevant ones are Umklapp pro- 
cesses. They correspond to backscattering of pairs of 
electrons and lead to the decay rate that scales as 7„ ^ 
rp2K^+-i 2ero doping and becomes exponentially small, 
~ exTp{—2kFVF/T), at large doping, see Eq. (27 1. This 



rate can be expressed in terms of the short distance cutoff 



^33^ and Umklapp coupling constants, see Eq. ( 27 ) or in terms 



of the cutoff-independent Umklapp gap, as in Eq. (31). 

In contrast to the Umklapp processes, the phonon- 
assisted decay rate depends on the doping level only 
on a very large scale vpkp ^ Tvp/sr ~ 0.6 x 10^ T. 
With the exception of rather low temperatures the prac- 
tically achievable doping levels lie below this scale. In 
this regime the phonon-assisted decay rate scales as 7ep ~ 
rp(i+Ka+)/2 ^ as given by Eq. (34). It can also be expressed 
in terms of a cutoff-independent energy scale Agp asso- 



ciated with electron-phonon interactions, see Eq. (38 1 



This scale is given by the gap in the electron spectrum 
induced by the twist Peierls instability. 

If the plasmon frequency is smaller than the tempera- 
ture, the ratio of the phonon-assisted rate and the Umk- 



The 



phonon assisted decay is dominant in the wide range of 
temperatures, T > (A^/Agp)^/"^. With the exception 
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of the region of very low temperatures it is dominated 
by electron-phonon scattering and scales as t'^+^^+^Z^. 
This is consistent with the temperature dependence of 
phonon-induced resistivitjff^. 
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Appendix A: fusion rules 



Equation (16) can be derived by noting that both the 



Uniklapp and electron-phonon backscattering Hamilto- 
nian densities have exponential dependence of the charge 
field <i>c+i ~ (g*ai7*<:+ -\-h.c.) and using the following iden- 
tity that holds for Gaussian averages, 



-(m^)2(0,+ (x)0,+ (a;i))o(e^'^''^^+(-i)e-^«''^^+(-^))o(0c+(2;2)0c+(2:'))o 
+(za^)2(0,+ (x)</)c+(xi))o(e»''*=+(-i)e-^'^"'^=+(-^))o(</)c+(xi).^c+(a:')>o + 

same terms with xi ^ X2 + ((/>c+(x)0c+(a;')>o(e"''''^'=+^'''^e-'"'''^<=+(^2))o. (Al) 



This formula can be established by expanding the ex- 
ponentials in the left hand side in the Taylor series and 
applying Wick's theorem to each term. The resulting 
series sums to the expression in the right hand side. 



The last term in the right hand side corresponds to 
a disconnected diagram and does not contribute to the 
correlation function. 
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